Approach to equilibrium of diffusion in a logarithmic potential 
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The late-time distribution function P(x,t) of a particle diffusing in a one-dimensional logarithmic 
potential is calculated for arbitrary initial conditions. We find a scaling solution with three surprising 
features: (i) the solution is given by two distinct scaling forms, corresponding to a diffusive (a; ~ t 1 / 2 ) 
and a subdiffusive (a; ~ V with a given 7 < 1/2) length scale, respectively, (ii) the overall scaling 
function is selected by the initial condition, and (iii) depending on the tail of the initial condition, 
the scaling exponent which characterizes the scaling function is found to exhibit a transition from 
a continuously varying to a fixed value. 

PACS numbers: 05.40.-a,05.10.Gg 



X 



I. INTRODUCTION 

There exist many physical systems whose temporal 
evolution is described by a diffusion process in a one- 
dimensional logarithmic potential. Examples include the 
denaturation process of DNA molecules [lj and the tem- 
poral evolution of the momentum distribution of cold 
atoms trapped in an optical lattice Effective log- 

arithmic potentials also appear, for example, in models 
of real-space condensation such as the zero-range process 
, relaxation to equilibrium of long-range interacting 
gases @, [9| , dynamics of the two-dimensional XY model 
below the Kostcrlitz-Thoulcss transition [l0[, the ABBM 
model for Barkhausen noise flU ] , and d yna mics of sleep- 
wake transitions during a night's sleep fl2j . 

The Fokker-Planck equation corresponding to diffusion 
in a one-dimensional potential is 
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where P(x, t) is the probability distribution. We consider 
potentials which increase logarithmically at large x, 



V(\x\ » 1) ~ 6Iog(|x|), 



(2) 



and are regular at x — 0. For simplicity, we have taken 
in these equations the diffusion constant and the temper- 
ature to be equal to 1. Here we consider the case b > 1, 
for which the system evolves into a stationary state given 
by the normalizable Boltzmann distribution 
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where Z = J e~ v ^d x is the normalization constant. For 
some applications the variable x is by definition non- 
negative, x > 0, as in the case of DNA denaturation 
where x corresponds to the length of a denaturated loop. 
In these cases, the equation has to be supplemented by 
a boundary condition at the origin. 

In this paper, we use a scaling analysis to study the 
long-time evolution of the probability distribution to- 
wards the stationary state. We find that the solution 



of Eq. ([1} relaxes to equilibrium via a universal scaling 
form which depends on the potential only through its 
asymptotic form ((2|). This scaling form exhibits several 
features which are not typically found in scaling solu- 
tions ff^j. (i) At large times, the equation exhibits two 
distinct scaling regimes which we refer to as the large- x 
and the small- a: regimes. The two scaling functions yield, 
to leading order in time, the distribution at any point x. 
They join smoothly at an intermediate scale xi(t), which 
grows with time, (ii) The overall scaling solution (com- 
posed of both regimes) is not unique. There exists a one 
parameter family of such solutions and the appropriate 
solution is selected by the tail of the initial distribution. 
The mechanism by which the initial condition selects the 
eventual scaling solution is analogous to that encountered 
in, e.g., fronts propagating into unstable states [H, [Hj]. 
(iii) For a class of initial conditions whose tails are suffi- 
ciently close to the eventual steady-state distribution, the 
scaling solution is found to be independent of the details 
of the initial condition. On the other hand, the scal- 
ing function resulting from other initial conditions varies 
continuously with the initial condition. 

A large- a; scaling solution of Eqs. {I])-© has recently 
been analyzed in 0, EH , where the dependence on the 
initial distribution has not been considered. Although 
this analysis is valid for a rather broad class of initial con- 
ditions, including compactly supported ones, other initial 
conditions which may arise in various physical circum- 
stances are left out. The analysis presented here applies 
to all initial conditions, and provides the scaling form in 
both the small- a; and large- a; regimes. 

The paper is organized as follows. In Sec. [HI wc present 
our scaling analysis. The scaling solution for the case of 
reflecting boundary conditions at the origin is considered 
in Sec. Ill A[ where a one-parameter family of solutions is 
identified. In Sec. Ill Bl we present a mechanism by which 
a particular solution is selected by the initial condition, 
and test it numerically. A generalization to other bound- 
ary conditions is discussed in Sec. Ill CI Finally, in Sec. 
IIIII we summarize our results. 
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II. SCALING ANALYSIS 



Scaling solutions for reflecting boundary 
conditions 



We begin our analysis by introducing a function G(x, t) 
defined as 



P(x,t)=P*(x)[l + G(x,t)}. 



(4) 



The deviation P*(x)G(x, t) from the steady state also 
satisfies Eq. (!]), but its normalization is zero. For x ^> 1, 
where the deviation of the potential from the logarithmic 
form © is negligible, G satisfies 



dG(x,t) bdG(x,t) d 2 G(x,t) 



dt 



dx 



dx 2 



(5) 



To be specific, we consider the case where the equa- 
tion is defined on the positive real axis, x > 0, with 
a reflecting boundary condition at x = 0. This implies 
8G{x = 0, t)/dx = 0. We later comment on other bound- 
ary conditions. 

The scaling solution of Eq. (5j) may in fact be ob- 
tained by an exact solution of Eqs. ([I])-©. This is 
done by transforming the Fokker-Planck equation into 
an imaginary-time Schrodinger equation via the trans- 
formation P(x,t) = e- vi - x ^ 2 ip(x,t) [jj]. The resulting 
equation for the "wavefunction" tp is 



dip(x,t) d 2 ip{x,t) 



dt 



dx 2 



V a (x)ip(x,t) 



(6) 



with the Schrodinger potential V 8 (x) = V ^ V ^ ■ 

For a potential of the form © this gives 



V s (x^> 1) ~ 7 /x 2 



(7) 



with 7 =5(5 + 1)- For large x, this equation describes 
the well studied problem of a quantum particle moving 
in a repulsive inverse square potential [17] . The solution 
of this equation may be found by expanding it in eigen- 
functions, ip(x,t) = J dkake~ k t ipk(x), where ?/>fc(x) can 
be expressed in terms of Bessel functions. The long-time 
scaling behavior is then obtained by studying the small-fc 
behavior of the amplitudes Ofe. Carrying out this expan- 
sion, we find that any localized initial condition for the 
function G evolves at long times to 



/ x 



6+i e -« 2 /4 



G{x,t) ~ T 1 f x [^=} with h{u)=u 

(8) 

The analysis is rather lengthy and will be presented else- 
where [l8(. Below we derive the long-time solution di- 
rectly by assuming a scaling form. We make use of the 
result (5J) only in relating the appropriate scaling solution 
to the initial condition. 

Let us now present the scaling solution of this equation 
and outline its derivation. As we demonstrate below, 
two length scales emerge at large times: a large- a; regime 



with x ~ i 1 / 2 , and a small- x regime, x ~ f , with a in- 
dependent 7 satisfying 7 < 1/2. Starting with the small- 
x regime, we consider a scaling solution of the form 



G(x,t)=t- 6 g(z), z=- 



(9) 



with some function g and exponents 7 and 5. Substitut- 
ing © in Eq. (5J) yields 

g"(z) - -g'(z) = -hzg'(z) + Sg{z)]t-^\ (10) 
z 

The right-hand side of this equation may be neglected as 
long as 7 < 1/2, yielding the solution 



g{z) = C + Cz b+ \ 



(11) 



where C and C are integration constants. Thus, for z <C 
1 one has G(zt 7 ,t) ~ t~ s C. In fact, since for z <C 1 this 
solution satisfies the boundary condition at x = 0, it is 
valid down to x = 0. 

In order to determine 7 and 5 one needs the solution 
in both small- a; and large- a; regimes. We thus consider a 
different scaling function for x ~ t 1 / 2 , 



G(x,t)=t-Pf(u), u = -^ 



x 



(12) 



where the scaling exponent j3 and the function f(u) are 
to be determined. Substituting (fT2"j) in Eq. ([5]) yields a 
family of ordinary differential equations for f(u), param- 
eterized by /3: 



en , ( u b \ ti 



/"+ U- + = o- 



(13) 



Requiring that the small- a; and large- a; solutions join 
smoothly at an intermediate scale 

F«a; 1 (t)«t 1/2 , (14) 

one concludes from (j9)), (fTTj) and (TT21 that 

/(«« 1) ~ Cu b+ \ (15) 

and 

£ = /3+(&+l)(i-7). (16) 
The solution of Eq. (JT3J) which satisfies (T5J is [l9[ 

t( \ r b+i P f l + b + 2P 6 + 3 u 2 \ 

f(u) = Cu + 1F1 I ;— ^— ;— 4)) ( 17 ) 

where ii*x is the hypcrgcometric function. For small and 
large arguments f(u) satisfies [lj| 

i b+1 for u < 1 

i)^{Bu~ 20 for u > 1,0 + 1 1,2,3,... , 

Bu^-^-t for u> l,/3 = 1,2,3,... 

(18) 



3 





FIG. 1. (Color online) A schematic representation of the 
scaling solution G(x,t) (see Eq. (|21l) ~) at a given late time 
( > 1 (not drawn to scale). The red double line represents the 
small- a; scaling form g(x /t 1 ^ b+1 ' l )t~ 6 while the blue solid line 
represents the large- a; scaling form f{x/t 1 ^ 2 )t~ 13 . The interval 
on which the two solutions overlap (Eq. (0) is labeled xi(t). 



FIG. 2. (Color online) The (a) small- a; and (b) large- a; scaling 
collapse obtained by numerically integrating Eq. (|27[) for b = 
3.5 and a = 1. Different curves correspond to different times 
in increasing order in the arrow direction. The dashed curve is 
the theoretical function (|21[) . with no fitting parameter. Inset 
shows a magnified region around the origin. 



where B is a known constant which depends on b and /?. 

Once the small- a; and large- x scaling functions are 
known, conservation of probability enables us to de- 
termine the scaling exponent 7 and the integration 
constant C. We evaluate the normalization condition 
J P*(x)G(x, t)dx = by splitting the integral into two 
domains, < x < xi(t) and x\(t) < x < 00. Using the 
solutions found in the two domains and evaluating the 
integrals to leading order in t we obtain 



1 



Also, 



7 



G 



6+1 



and 6 = (3 + 



1 



(19) 



c = -zJ ' u - b ^ du = 4-W^T- (20) 

Summarizing the results of the scaling analysis, we find 
that the solution of Eqs. (H}-© is, to leading order in t, 



P( X ,t) Si P*(x) +P*(x) 



—5 



X < X\(t) 



f(&)t-f>, x> Xl (t), 



(21) 

where / and g are given in Eqs. (fTTj) and ([T7|). 7 and S 
are given in (fH?)) . and Xi(t) satisfies (fl"4")l . We thus obtain 
a one-parameter family of scaling solutions parameter- 
ized by /?. Below we denote a member of this family by 
Pfs(x, t) = P*(x)[l + Gfj{x, t)}. The overall form of these 
scaling solutions for G is given schematically in Fig. [T] 
For small-x, the function is flat, with a value that ap- 
proaches zero as t~ s . For x > x\(t). it exhibits a peak 
at x ~ s/t, whose height scales as t~@ . Thus, as time 
progresses, the peak shrinks and moves to the right, and 
G approaches its steady-state value G(x,t — > 00) = 0. 



B. Selection mechanism 

To conclude our analysis, we argue now that the pa- 
rameter j3 is selected by the tail of the initial distribution. 



To this end, we consider initial distributions of the form 



G (a; > 1) ~ Ax' 



(22) 



where a > b — 1 and A are parameters. Initial conditions 
Gq(x) which decay faster than algebraically correspond 
to a = 00. Note, however, that localized initial distri- 
butions P(x, 0) correspond to a = 0. Since the dynami- 
cal process of Eq. ([1]) is diffusive, one would naively ex- 
pect the tail of the distribution to remain unchanged for 
x 7S> \fi. Using the asymptotics (fT8|) . this suggests that 
f3 = a/2. The parameter C is then given by C = A/B. 
Substituting the value of B yields [19| 



rq-q/2) 

Ol+fr+ap( 3+fc ^ 



(23) 



This naive argument is found to be valid only for a < 2. 
For a > 2 we make use of Eq. |(5J) to demonstrate that the 
correct scaling solution is given by (3 = 1. This is done by 
analyzing the stability of a scaling solution with a given 
/J. We thus consider a perturbation SP(x,t) around the 
scaling solution ([2~Tj) : 



P(x,t)-P*(x) =P*(x)Gp(x,t)+5P(x,t), (24) 

which is initially localized in x, such as a function with 
compact support. Normalization dictates that J SPdx = 
0. This perturbation satisfies Eqs. (HJ)-©. 

The scaling solution Gp is stable to such perturbations 
and will dominate the approach to equilibrium only if at 
late times SP(x, t) is negligible compared to it. The exact 
solution of Eq. (QJ reveals that the scaling solution © 
to which a localized initial condition converges at long 
times is of the form (|2"Tj) with /3 = 1 [lj|. This can be 
understood heuristically by noting that /3 = 1 is the most 
localized of all scaling solutions, see (fT5]l . Therefore, the 
scaling form at large- a; of Eq. (f2"4"| is given by 



(25) 



P(x, t) - P* (x) = P* {x)G (x,t) + 5P{x, t) 



P*(x) 



1/1 (71) 
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where fp is the solution (|T7)) corresponding to f3. For 
f3 < 1 the second term on the right-hand side is negligible 
compared to the first, and therefore the fp solution is 
stable. On the other hand, for j3 > 1 the second term is 
dominant and the scaling solution is given by f\. Thus, 

The exact solution of the Fokker-Planck equation further 
reveals that when a = 2 there are logarithmic corrections 
to Eq. d2TJ [H[. Note that for a > 2, the constant C is 
not given by (|23[) and it depends on the details of the ini- 
tial condition. The large- a; scaling solution for localized 
initial distributions (corresponding to a = 0) agrees with 
the results previously obtained in 0, 0, HH • Knowledge 
of the scaling solution for other initial conditions is often 
required for calculating correlation functions of physical 
interest. Applications of this approach to physical exam- 
ples where the late time behavior is determined by the 
initial conditions will be presented elsewhere fl8| . 

In order to check the applicability of the scaling solu- 
tion found in this analysis, we studied numerically the 
evolution of a single-site zero-range process at critical- 
ity In this process, particles hop into a site with 
a constant rate 1, and hop out of this site with rate 
w(n) = 1 + b/n, where n is the number of particles in 
the site. The occupation number probability distribu- 
tion P(n, t) satisfies the master equation 

Q 

— P(n) = P(n-l) + w(n+l)P(n+l) - [1 + w(n)]P(n). 

(27) 

This is a discrete version of (JI)-©. Its steady state is 
given by P*(n) oc [w(l) . . . iu(n,)] -1 ~ n~ b . We have 
studied the relaxation to the steady state starting from 
various initial conditions. In Fig. [2] we display the results 
obtained for an initial condition with a = 1. A very 
good agreement with the predicted scaling functions is 
obtained at long times. To investigate the dependence 
(|2T))) of the scaling exponents on the initial condition, we 
have measured G(0, t) which is predicted to decay to zero 
as t~ s (fl9|) . As an independent measure, we calculated 
(n(t)), which is predicted to decay to its equilibrium value 
(n)eq as t~ a with a = (3 + 6/2 — 1. This can be easily 
obtained using the scaling function (|2T|) : 

/>oo 

(n(t)> - (n)eq « t-^+ 6 /a-i)c / u l - b f{u)du. (28) 

Jo 

In Fig. [3] we compare numerical measurements of 6 and a 
to the theoretical predictions and find a very good agree- 
ment, both for a < 2 and a > 2. 

C. Other boundary conditions 

The analysis presented above may be extended to con- 
sider other boundary conditions as well. For example, 
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FIG. 3. The scaling exponents 5 and a obtained by numeri- 
cal integration of (|27[) , compared with theoretically predicted 
values. Here b — 2.9 and A = 1. The transition at a — 2 is 
readily seen. 

for an absorbing boundary condition at x = 0, namely 
P(0,t) = 0, the stationary distribution vanishes. Defin- 
ing G(x, t) by 

P(x,t) =e- v{x) G{ Xl t) (29) 

and repeating the derivation outlined above for the evo- 
lution of G, we find that in the scaling limit it takes the 
same form (|12[) as before, with / unchanged, 7 = and a 
different function g(z) which can be calculated explicitly 
[HI ■ Equation (f2l)]) , which yields the relation between ft 
and a, remains valid. However, since G is defined dif- 
ferently in this case, the same initial distribution P(x, 0) 
would correspond to different values of a depending on 
whether the boundary condition is reflecting or absorb- 
ing. Other boundary conditions, and the case where the 
equation is defined on the entire x-axis, may be treated 
similarly flU ]. 

III. CONCLUSION 

In summary, the approach to equilibrium of a diffusion 
process in a logarithmic potential is analyzed in the scal- 
ing limit and is shown to exhibit uncommon and inter- 
esting features. These include the existence of two char- 
acteristic length scales, the fact that the scaling function 
depends on the tail of the initial distribution, and the 
non-analytic way the scaling exponents varies with the 
initial condition. The mechanism by which the scaling 
exponents are selected is similar to the one encountered 
in p roblems of velocity selection in propagating fronts 
[l3l [Til]. Here, however, unlike problems of front prop- 
agation, the evolution equation is linear, although inho- 
mogeneous. This facilitates the explicit demonstration of 
the selection mechanism. 
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